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Abstract 

The rare decay B + — > Dfcp can occur only via annihilation type diagram 
in the standard model. The small branching ratio predicted in the standard 
model makes this channel sensitive to new physics contributions. We analyze 
this decay mode, both in the standard model and in several extensions of 
it. The models considered are minimal supersymmetric model with R-parity 
violation and two Higgs doublet model. The experimental verification of our 
findings of large branching ratio and/or nonzero CP asymmetry may signal 
the presence of new physics. 
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Recent results from the ongoing experiments in 5-physics at BABAR and BELLE have 
attracted a lot of attention. The main objective of these B experiments is to explore in detail 
the physics of CP violation, to determine many of the flavor parameters of the standard 
model (SM) at high precision and to probe for possible effects of new physics beyond the 
SM [T|-0|. The intensive search for Physics beyond the standard model is performed now a 
days in various areas of particle physics. The B system thus offers a complementary probe 
to the search for new physics. In B experiments, new physics beyond the SM may manifest 
itself in the following two ways : 

i. decays which are expected to be rare in the SM and are found to have large branching 
ratios. 

ii. CP violating asymmetries which are expected to vanish or to be very small in the SM 
are found to be significantly large. 

Thus the rare B meson decays are suggested to give good opportunities for discovering 
new physics beyond the SM. Since their branching ratios are small in the SM, they are very 
sensitive to new physics contributions. In the last few years, different experimental groups 
have been accumulating plenty of data for the rare B decay modes. Some of them have 
already been measured by the B factories in KEK and SLAC. 

In this context it is interesting to analyze the rare decay B + — > D+(f), which is a pure 
annihilation type decay in the SM. The four valence quarks in the final states D s and </> 
are different from the ones in the parent B meson, i.e. there is no spectator quark in this 
decay. In the usual factorization approach, this decay mode can be described as the b and 
u quarks in the initial B meson annihilating into vacuum and the final D s and (f) mesons 
are produced from the vacuum afterwards. The dynamics of exclusive hadronic B decays 
occurring via the W- exchange or annihilation diagrams, is not yet understood. The decay 
rates for such transitions are argued to be negligibly small due to the suppression of helicity 
and (or) form factors. However, a solid justification of this argument is necessary in both 
theory and experiments. This decay mode has been recently studied in perturbative QCD 
approach with a branching ratio 3 x 10~ 7 @], which is far below the current experimental 
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upper limit 

Br(B + -> D+(j)) < 3.2 x 10~ 4 (1) 

Therefore it provides an appropriate testing ground for physics beyond the SM. 

In this paper we intend to study the decay mode B + — > Df<p, both in the standard model 
and in several extensions of it. The models considered here are the Two Higgs Doublet Model 
(2HDM) |J and minimal supersymmetric model with R-parity violation (RPV) 0. 

We first consider the contributions arising from the SM, where the effective Hamiltonian 
describing the decay mode is given as 

n eff = ^v: b v cs [dMOM + cmom] (2) 

where 

1 = [6 7 m(1-7 5 )H[c7 m (1-7 5 )s] 

2 = [h,(l -T5)s][5T M (l -7sH (3) 

and 6*1,2 (//) are Wilson coefficients evaluated at the renormalization scale fi. In the next to 
leading logarithmic approximation their values are evaluated at the b quark mass scale as 

§ 

C x {m b ) = 1.082 and C 2 (m b ) = -0.185 (4) 

The corresponding transition amplitude is given as 

Amp( J B+ - D S <P) = ^(D^Kc^il ~ J 5 )s){h»(l ~ 1»\B + ) (5) 

where a± — C\ + C^/A^,, with iV<7 is the number of colors. The evaluation of the matrix 
element of four fermion operator from the first principles of QCD is an extremely demanding 
challenge. To have some idea of the magnitudes of matrix element, one can still use the 
factorization method, factorizing the four quark operators relevant to non leptonic B decays 
into the product of two currents and evaluating separately the matrix elements of the two 



currents. Thus in the factorization approximation one can write the corresponding transition 
amplitude as 



Amp(B+-Z>+0 = -^V: b V cs a 1 <Z>+0|(&y*(l- j 5 )s) \0}(0\(h,(l ~ l 5 )u)\B + ) (6) 



The main uncertainties in evaluating the transition amplitude are due to the matrix element 
like (.0+01(07^(1 — 7 5 )s)|0) with (p Dg + p<f,) 2 — M§. To get an idea about the contributions 
from annihilation diagrams, one can assume single pole dominance for the matrix element 
and relate it to the crossed channel (</>|s7 M (l — j5)c\D s ). We see that the annihilations do 
not give a large contributions here : with single pole dominance, form factors are suppressed 
by a factor M D 2 J(M B — M^) because the transferred momentum ^[Q 1 = Mb is large with 
respect to pole mass . The annihilation contribution could be larger if there were other 
pseudoscalar (cs) resonators heavier than the pole mass Mp , which would enhance the form 
factors. The matrix element of the pseudoscalar and vector meson is usually decomposed as 



where Q = (p — q) is the momentum transferred during the transition process. Using the 
decay constant relation 



i 



(<f>( q ,e)\s^(l- l5 )c\D s (p)) 




-(p + q yA 2 {Q 2 ) 



(7) 



(0\b^j 5 u\B 



(Q)) = if B Q 



(8) 



the transition amplitude in the SM is given as 



Amp(5 -> D s ct)) 



V2 



V: b V cs ai/fl 2M 4> A {M B ){e* ■ Q) 



(9) 



and the corresponding decay width as 




\Am P (B^D s <P)/(e*-Q)\ 2 , 



(10) 
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where p c is the cm. momentum of the decay particles. The value of the form factor Ao(0) 
at zero momentum transfer for the transition D s — > <fi is given as || v4 (0) = 0.7 and the Q 2 
dependence of the form factor is given as 

MQ 2 ) = —V ( n ) 

where Mp is the pole mass with value Mp = 1.97 GeV |§. Thus we get 

Aq(M 2 b ) = -0.113 (12) 



Using ai=1.04 [10[, which is extracted from the experimental data on B — ► D*(ir, p), /#=190 
MeV, the CKM matrix elements from || and Tp = 1.653 x 10~ 12 s, we obtain the branching 
ratio in the SM as 

Br(B + -> = 1.88 x 10~ 6 (13) 

Thus one can see that the obtained branching ratio in the SM using the factorization as- 
sumption is far below the experimental value Br(B + — > D+cf)) < 3.2 x 10~ 4 . Furthermore, 
it should be noted here that the direct CP asymmetry for this decay mode is zero in the SM 
since it receives contribution only from the single annihilation diagram. 

We now proceed to evaluate the branching ratio in the QCD Improved factorization 
method, which has been developed recently to study the hadronic B decays. This 



method incorporates elements of naive factorization approach as its leading term and per- 
turbative QCD corrections as subleading contributions and thus allowing one to compute 
systematic radiative corrections to the naive factorization for hadronic B decays. This 
method is expected to give a good estimate of the magnitudes of the hadronic matrix el- 
ements in non leptonic B decays. However, in the QCD factorization approach the weak 
annihilation contributions are power suppressed as h.QCD/ m b and hence do not appear in the 
factorization formula. Besides power suppression they also exhibit end point singularities 
even at twist two order in the light cone expansion of the final state mesons and therefore can 
not be computed self consistently in the context of hard scattering approach. One possible 



way to go around the problem is to treat the end point divergence arising from different 



sources as different phenomenological parameters [T2|]. The corresponding price one has to 
pay is the introduction of model dependence and extra numerical uncertainties. In this 
work we will follow the treatment of |12j and express the weak annihilation amplitude for 



B + -> D+6 as 



Amp ( J B+ - Dt<f>) = ^V: b V cs f B f D JM D s, 0) (14) 
The annihilation parameter bi(D s ,<p) is given as 

h{D.,<f>) = jZc 1 A[{D.,<l>) (15) 

where the color factor Cp = ^ 1 and N c = 3. The function A\(D S , 0) is given as 

fi T i i 

+ ± , (16) 



»1 rl 

A\{D S ,4>) = ira s I dx dytp Ds (x)^(y) 



y(l - xy) x 2 y 

where ipD s {x) and ip<p(y) are the light cone distribution amplitudes (LCDA) for the final 
mesons and x is the longitudinal momentum fraction of c quark in D s and y is the momentum 
fraction of s in 0. Using the assumption that the LCDAs of the mesons D s and <fi are 
symmetric, one can parameterize the weak annihilation contribution as [T3| 

A\(D S , 0) ~ 187ra s (X A + y - 4) (17) 

where Xa = Jq dx/x parameterizes the end point divergence as 

f 1 dx , M B ,a 
X A = / — = \n^- + pe- ie (18) 
Jo x A 

p varies from to 6 and 9 is an arbitrary phase < 9 < 360°. Using A = Aqcd = 200 
MeV and pe id = m as default values along with f B = 190 MeV, / = 233 MeV, f Ds = 280 
MeV, Ci(mi,) = 1.082 and a s {mi) = 0.221, we obtain the branching ratio in QCD improved 
factorization approach as 

Br(B + -> D+(f>) = 0.67 x 10" 6 (19) 
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We now proceed to calculate the branching ratio for this decay mode in type II of Two 
Higgs Doublet Model 0. In type II 2HDM, the up-type quarks get mass from one doublet, 
while down- type quarks and charged leptons from the other doublet. Charged Higgs Yukawa 
couplings are controlled by the parameter tan/5 = v^/vi, the ratio of vacuum expectation 
values of the two doublets, normally expected to be of order m t /m&. For our concern the 
if ± effectively induce the four fermion interaction as 



nj2HDM 



Gf Vub V cs I - 



m h X{\ + 75) + m u Y{l - 75) 



X < s 



m c F(l + 7 5 )+m s X(l-7 5 ) 



(20) 



where m# denotes the mass of the lightest charged scalar particle, m g 's denote the con- 
stituent quark masses and X ~ 1/Y = 1)2/ V\ = tan/3. In the above equation the terms 
proportional to Y can be safely neglected as X is generally taken as large flUf . To evaluate 
the matrix elements we use the equation of motion to transform the (S — P)(S + P) currents 
to corresponding (V — A)(V + A) form as 



±75)92 



^(gi7/xg2) ^(gi7 M 7 5 g 2 ) 



(21) 



m q2 + m qi 

where the quark masses are current quark masses. Thus we obtain the expressions for matrix 
elements as 



(0\(Hu)\B + ) 



if B M_ 



B 



(22) 



and 



\sj 5 c\D s ) = i(e ■ Q)- 



2M A 



-A Q {Q 2 ) . 



m s + m c 

Using the above relations, one can obtain the transition amplitude in 2HDM as 



(23) 



Amp(£+ - DM*hdm = %V^vJ^f B 



mbfn s M B 



V2 



M 



H 



Q)2M (j) A (Q 2 ) (24) 



So the total amplitude in 2HDM including SM contributions 



A SM+2HDM = A SM^ _ Ri] 



(25) 



where 

1 tan 2 /3 m b m s M B 

ti,\ — —75 — ; (^b) 

ai M H (m b + m u ){m s + m c ) 

The free parameters of the 2HDM namely tan/3, and Mh are not arbitrary, but there are 
some semi quantitative restrictions on them using the existing experimental data. The most 
direct bound on charged Higgs boson mass comes from the top quark decays, which yield 



the bound M# > 147GeV for large tan/? [yj|. Furthermore, there are no experimental 



upper bounds on the mass of the charged Higgs boson, but one generally expects to have 
Mh < ITeV in order that perturbation theory remains valid [|16|]. For large tan/5 the most 



stringent constraint on tan/5 and Mh is actually on their ratio, tan/3/M#. The current 
limits come from the measured branching ratio for the inclusive decay B — > Xri>, giving 
tan P/Mh < QAQGeV^ 1 [fL4| . Using the constituent quark masses as ra& = 4.88 GeV, 



m s = 0.5 GeV [[T7[] and current quark masses at the b quark mass scale as m b = 4.34 GeV, 
m c = 0.95 GeV, m s = 90 MeV and m u = 3.2 MeV Q, the branching ratio in 2HDM is 
found to be 

Br {B + -> Dt<j))\2HDM = 8.0 x lO" 6 (27) 

We now analyze the possibility of observing direct CP violation in this decay mode 
since it now receives contributions both from the SM and 2HDM. We can write the decay 
amplitude ( p5|) as 

A SM+2HDM = \ A SM^ /+ 8 l} _ \ Rl ^(<P+^ ^ (2g) 

where 7 = arg(V* b ) and 8\ are the weak and strong phases of standard model amplitudes. 
4> and S are the relative weak and strong phases between 2HDM and SM amplitudes. Thus 
the direct CP asymmetry for the decay mode is given as 

Br(B+ -> D+cf)) - Br(B~ -> L> 6 70) 

a, — 



cp Br(B+ -> D+(p) + Br(B- -> D^^) 

(29) 



2|i2i 


sin 


sin 5 


1 + 




2 _ 


2\Ri 


COS cos 5 



If we set = 5 = 7r/2, the maximum possible value of direct CP violation in 2HDM is found 
to be 



a cp (B 




<t>)\2HDM <59% 



(30) 



We now analyze the decay mode in minimal supersymmetric model with R-parity violation. 
In the supersymmetric models there may be interactions which violate the baryon number B 
and the lepton number L generically. The simultaneous presence of both L and B number 
violating operators induce rapid proton decay which may contradict strict experimental 
bound. In order to keep the proton lifetime within experimental limit, one needs to impose 
additional symmetry beyond the SM gauge symmetry to force the unwanted baryon and 
lepton number violating interactions to vanish. In most cases this has been done by imposing 
a discrete symmetry called R-parity defined as R p = (— l)( 3 - B+L + 25 ) ) where S is the intrinsic 
spin of the particles. Thus the -R-parity can be used to distinguish the particle (R p =+1) 
from its superpartner (R p = — 1). This symmetry not only forbids rapid proton decay, it also 
render stable the lightest supersymmetric particle (LSP). However, this symmetry is ad hoc 
in nature. There is no theoretical arguments in support of this discrete symmetry. Hence it 
is interesting to see the phenomenological consequences of the breaking of R-parity in such 
a way that either B and L number is violated, both are not simultaneously violated, thus 
avoiding rapid proton decays. Extensive studies has been done to look for the direct as well 
as indirect evidence of R-parity violation from different processes and to put constraints 
on various R-parity violating couplings. The most general -R-parity and Lepton number 
violating super-potential is given as 



where, i,j, k are generation indices, Li and Qj are SU(2) doublet lepton and quark su- 
perfields and E%, D% are lepton and down type quark singlet superfields. Further, is 
antisymmetric under the interchange of the first two generation indices. Thus the relevant 
four fermion interaction induced by the R-parity and lepton number violating model is 



W V = i^ijkLiLjEl + X^LiQjDl , 



(31) 
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n * = " (1 " 75)6 s_(1 + 7s)c (32) 

where the summation over i = 1, 2, 3 is implied. It should be noted that the RPV Hamilto- 
nian has the same form as the 2HDM Hamiltonian except the couplings. So using the Eqs. 



3) and ( |23|) one can easily obtain the amplitude for the B + — > Dfcf) in RPV model 

4:M~ L (m b + m u ) (m s + m 



Amp(5+ - DtMsry = ■ }b ™ b ■ ,(e • Q)2M,A (Q 2 ) (33) 



and the total amplitude as 

^SM+W = _ ^ ^ (34) 

where 



R 2 = - 7; 2 f2 il3 /m ; Z ; Z , (35) 



v 7 ^ 1 A^ASa 

G^KTb^ai 4M| (m& + m u )(m s + m. 

Using A 2i2 A^ 3 = 2.88 x 10" 3 [0, we obtain the branching ratio in RPV model as 



Br(B + -> £>+0) = 3.06 x 10~ 4 (36) 

Thus we found that the branching ratio in RPV model is quite large in comparison to SM 
prediction. So if the experimental value is found to be in this range, it will definitely be a 
signal of new physics beyond the SM. 

Replacing R 2 in place of R\ in Eq. fl2"9|) , the direct CP asymmetry in the decay mode 
B + — > in RPV model is found to be 

a cp (B + ^D+<J))\ R p V <U% (37) 

To conclude, in this paper we have calculated the branching ratio of the two body 
hadronic decay mode B + — > Df<f) in the standard model as well as in the two Higgs doublet 
model and RPV model. We found that the branching ratio in the RPV model is quite large 
in comparison to the SM prediction, whereas the 2HDM prediction is approximately one 
order higher than the SM value. The direct CP asymmetry a cp , which is expected to be zero 
in the SM, is found to be nonzero and large in both RPV and 2HDM analyses. From our 
analyses it follows that the new physics contribution is quite large and significant. Therefore 
the rare decay mode B + — > D^<p provides an ideal testing ground to look for new physics. 
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